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Abstract 

Background: Maize rough dwarf disease (MRDD) is a devastating viral disease tliat results in considerable yield 
losses worldwide. Three major strains of virus cause MRDD, including maize rough dwarf virus in Europe, Mai de Rio 
Cuarto virus in South America, and rice black-streaked dwarf virus in East Asia. These viral pathogens belong to the 
genus fijivirus in the family Reoviridae. Resistance against MRDD is a complex trait that involves a number of 
quantitative trait loci (QTL). The primary approach used to minimize yield losses from these viruses is to breed and 
deploy resistant maize hybrids. 

Results: Of the 50 heterogeneous inbred families (HIFs), 24 showed consistent responses to MRDD across different 
years and locations, in which 9 were resistant and 15 were susceptible. We performed trait-marker association 
analysis on the 24 HIFs and found six chromosomal regions which were putatively associated with MRDD 
resistance. We then conducted QTL analysis and detected a major resistance QTL, qMrddl, on chromosome 8. By 
applying recombinant-derived progeny testing to self-pollinated backcrossed families, we fine-mapped the qMrddl 
locus into a 1.2-Mb region flanked by markers Ml 03-4 and Ml 05-3. The qMrddl locus acted in a recessive manner 
to reduce the disease-severity index (DSI) by 24.2-39.3%. The genetic effect of qMrddl was validated using another 
Eg recombinant inbred line (RIL) population in which MRDD resistance was segregating and two genotypes at the 
qMrddl locus differed significantly in DSI values. 

Conclusions: The qMrddl locus is a major resistance QTL, acting in a recessive manner to increase maize resistance 
to MRDD. We mapped qMrddl to a 1.2-Mb region, which will enable the introgression of qMrddl -based resistance 
into elite maize hybrids and reduce MRDD-related crop losses. 
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Background over three million hm^ of crops each year. Yield losses are 

Maize rough dwarf disease (MRDD) is a viral disease generally over 30% in affected areas and can reach 100% in 

that results in substantial yield losses in Europe, East regions of severe infection [6]. The virus that causes MRDD 

Asia, and South America [1-4]. MRDD was discovered belongs to the genus Fijivirus in the family Reoviridae, but 

in 1954 in China (South Xinjiang and West Gansu) and virus strains vary between continents. MRDD is caused by 

has posed a grave threat to maize production during the maize rough dwarf virus in Europe, Mai de Rio Cuarto virus 

last two decades, especially in the Yellow-Huai-Hai River in South America, and rice black-streaked dwarf virus in 

plain [5]. Between 2008 and 2011, MRDD has affected East Asia [7]. These viruses are transmitted in a persistent 

manner by planthopper insect vectors [8]. 
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controlling MRDD include pesticides, shifting the date(s) 
when seeds are planted (i.e., based on projected insect 
populations), and improving field management [9]. These 
methods limit the planthopper population and reduce, to 
some extent, MRDD severity, but always with high risk 
and low efficiency. The identification of MRDD -resistant 
strains, however, likely represents the most cost-effective 
and environmentally friendly way to minimize yield losses. 
It is therefore important to develop and deploy resistant 
hybrids by mapping and cloning genes and quantitative 
trait loci (QTLs) that confer resistance to MRDD. 

Under natural-infection conditions, the maize germ- 
plasm displays variable resistance to MRDD [10-14]. 
The major source of resistance is derived from US hy- 
brid P78599. Evaluation of 96 inbred lines and 136 hy- 
brids suggests that MRDD resistance is a quantitative 
trait [15]. Wang et al. (2000) reported that maize resist- 
ance to MRDD is a quantitative trait controlled by many 
genes, each with a small effect [16]. In Argentina, a par- 
tially resistant line yielded moderate heritability of resist- 
ance to the Mai de Rio Cuarto virus, ranging from 0.44 
to 0.56 [4]. Using an ¥23 QTL-mapping strategy, two 
QTLs were identified on bins 1.03 and 8.03/4 that to- 
gether explained 36.2% of the phenotypic variance [17]. 
A major QTL on chromosome 8 for MRDD resistance 
was identified in the Chinese maize inbred line, XI 78, 
based on 514 gene-derived single nucleotide polymor- 
phisms (SNPs) [6]. Using an F2 population derived from 
the highly resistant line 90110 and the susceptible line 
Ye478, Luan (2012) found at least three QTLs within 
chromosome bins 6.02, 7.02, and 8.07 that confer MRDD 
resistance [18]. 

In this study, we applied trait-marker association to 24 
heterogeneous inbred families (HIFs) and QTL analysis 
to segregating population derived from HIFs to identify 
regions of the maize genome that affect resistance to 
MRDD. We then fine-mapped the major QTL by subjecting 
self-pollinated backcrossed families to recombinant-derived 
progeny testing. Finally, an F^ recombinant inbred line (RIL) 
population was used to validate the effect of this QTL. 
These results provide valuable information concerning 
maize resistance to MRDD, and markers developed within 
the qMrddl region may prove useful in resistance breeding 
programs. 

Results 

Evaluation of HIFs In resistance to MRDD 

The 50 HIFs developed from a maize hybrid CL1165 
were evaluated for their resistance to MRDD in Taian 
for three years (from 2008 to 2010) and another two 
locations (Feicheng and Jining) in 2010. Of them, 9 
displayed consistent resistance regardless of the year 
or location; while 15 showed high susceptibility in 
every location across three years (Additional file 1). 



These 24 HIFs were therefore used for subsequent 
trait-marker analysis (Figure 1). 

Trait-marker association analysis In HIFs 

Genotyping 56,110 SNPs on the 24 HIFs generated 
48,384 successful calls based on a cluster-separation 
score of >0.3 and <50% missing data. Of these 48,384 
SNPs, 8,668 were polymorphic with a minor allelic fre- 
quency > 0.05 among the 24 HIFs. Importantly, 105 SNPs 
co-segregated with MRDD resistance. After projecting 
onto the maize accessioned golden path (AGPv2), 103 co- 
segregating SNPs were found to cluster within six 
chromosomal regions (Table 1). The remaining two SNPs 
likely resulted from random error and were excluded from 
further analysis. The six genomic regions represent candi- 
date MRDD-resistant loci. 

QTL analysis of maize resistance to MRDD 

Based on both genotypes and phenotypes, four HIFs, 
NT401 (susceptible), NT399 (resistant), NT409 (suscep- 
tible), and NT411 (resistant), were selected to prepare 
segregating populations. The Fi hybrid between NT401 
and NT399 was backcrossed twice to the susceptible 
NT401 to generate BC2F1 populations. The F^ hybrid 
between NT409 and NT411 was backcrossed to the sus- 
ceptible NT409 and then self-pollinated to generate 
BC1F2 populations (Figure 1). The totally 485 BC2F1 and 
211 BC1F2 families were then evaluated for MRDD sus- 
ceptibility at four experimental stations: Taian, Feicheng, 
Heze, and Jining. The disease-severity index (DSI) values 
were estimated for BCi individuals based on their BC2F1 
or BC1F2 families. The resultant DSI values displayed 
continuous distributions, ranging from 66.7-100% in 
Taian (BC1/BC2F1), 39.7-100% in Feicheng (BC1/BC2F1), 
51.9-100% in Taian (BC1/BC1F2), and 67.9-100% in 
Jining (BC1/BC1F2), respectively (Figure 2B), implying 
the quantitative nature of maize resistance to MRDD. 
Intriguingly, DSI distributions skewed toward the sus- 
ceptible parent, suggesting that a major dominant resist- 
ance QTL was absent. In Taian, for example, BC2F1 
individuals were extremely susceptible to MRDD, and 
BC1F2 individuals exhibited a segregation bias towards 
susceptibility (Figure 2B). In addition, location had a large 
effect on DSI. A population planted in Jining generally 
had a higher DSI than the same BC2F1 family in Taian. 
There was no MRDD outbreak in Heze, so nearly every 
family appeared MRDD resistant (even the susceptible 
line, NT409). Data from this location were therefore 
excluded from further analysis. As such, environmental 
conditions could fluctuate wildly in MRDD resistance 
(Figure 2B). 

Within the six candidate regions, 18 simple-sequence 
repeat (SSR) primer pairs were retrieved from the public 
database, and additional 81 primer pairs were developed. 
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Figure 1 Experimental flow chart for QTL identification and fine-mapping. Twenty-four HIFs were selected for trait-marker association 
analysis. Associated regions were validated using two BCi populations, and a major QTL, qMrddl, was identified and mapped. Recombinants 
identified within the BC1F4 population were genotyped and phenotyped to fine-map the QTL. 



These 99 SSR markers were used to screen for polymor- 
phisms between the parental lines, resulting in the iden- 
tification of 14 polymorphic SSR markers (Table 2 and 
Additional file 2). Each BCi individual was then geno- 
typed using these 14 SSR markers to distinguish hetero- 
zygous from homozygous genotypes at the six candidate 
regions. The phenotype of a given BCi individual was 



Table 1 Co-segregating MRDD resistance regions 
identified through trait-marker association 



Chr. 


Location 
Start point End point 


Number 
of SNPs 


Magnitude (bp) 


1 


33,678,550 


35,134,250 


21 


1,455,710 


3 


213,804,134 


214,476,941 


4 


672,807 


4 


216,519,340 


234,819,022 


48 


18,299,682 


5 


1 7,340,597 


24,441,850 


3 


7,101,253 


8 


114,594,605 


168,285,077 


10 


53,690,472 


9 


153,489,841 


153,784,681 


17 


294,840 



Location corresponds to maize B73 reference genome APG V2. 



represented by DSI values calculated from its BC1F2 or 
BC2F1 progeny. Mean DSI value was calculated individu- 
ally for either heterozygous or homozygous BCi geno- 
types at each of the six candidate regions. Based on DSI 
values from BC1F2 progeny at Taian and Jining, three 
SSR markers within bins 8.04/05, bnlgl62, umcl670, 
and umcll72, were significantly associated with MRDD 
resistance, as difference in DSI was significant between 
heterozygous and homozygous BCi genotypes (P < 0.01). 
Another two SSR markers within bin 5.03, C5-5 and 
C5-24, were marginally associated with MRDD resist- 
ance in Taian (P = 0.04 and P = 0.06, respectively) 
(Table 3). For bins 8.04/05, the mean difference in DSI 
between homozygous and heterozygous BCi genotypes 
differed by 9.6% and 7.0% in Taian and Jining, respect- 
ively. At the other five chromosomal regions, however, 
differences in DSI between these two BCi genotypes 
were less than 4.0%. Based on DSI values from BC2F1 
progeny at Taian and Feicheng, three SSR markers on 
bins 8.04/05 and two SSR markers on bin 5.03 were 
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Figure 2 Phenotype of MRDD. A. Plants with MRDD scores of 1-9. Disease scores were primarily based on plant height. Here 1 = no symptom; 
3 = slightly shorten superior internodes, about 4/5 plant height of healthy plant; 5 = dark-green leaves, waxy enations on abaxial surfaces of leaves 
and sheaths, obviously shorten superior internodes, about 2/3 plant height of healthy plant; 7 = severe shorten internodes, malformed tassels, 
about 1/2 plant height of healthy plant; 9 = severe stunning; suppressed flowering, and a lack of ears; plant height <l/3 of healthy plant. B. DSI 
distributions within segregating populations. Vertical arrows represent the DSI of parental lines. 



marginally associated with MRDD resistance (Table 3). 
The homozygous and heterozygous BCi genotypes 
exhibited very similar DSI values (<2.0% difference). 
Taken together, the bins 8.04/05 region was most likely 
to contain a major MRDD resistance QTL. 

The BC1F2 population was selected to map the resist- 
ance QTL since it showed stronger segregation in 
MRDD resistance than the BC2F1 population. Based on 
B73 reference sequence, 63 PCR-based markers were 



synthesized within bins 8.04/05 and nine of them were 
polymorphic between the parental lines, NT409 and 
NT411. Besides three previous SSR markers, these nine 
markers were also used to genotype the 211 BCi individ- 
uals. Based on genotypic data and the DSI values from 
BC1F2 families, we performed linkage mapping to define 
the resistance QTL (in Jining and Taian). A resistance 
QTL was detected within a 15-Mb region flanked by the 
markers M103-4 and umcll72 (Figure 3A). This resistance 
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Table 2 Markers developed to map the qMrddl locus 



Name 


Chr. 


Location (Mb) 


Forward primer 


Reverse primer 


Annealing 
temp, ro 


Type 


Ml 03-1 


8 


103.18 


AGAGGAGGCTOGATGCGGT 


TOAAAAAAAGGCATAGGCT 


60 


SSR 


Ml 03-4 


8 


103.93 


CACTCTCATCCTCTCGCACC 


ATCTAATCACAGCGCGCAAG 


60 


STS 


Ml 03-7* 


8 


104.17 


AACGAGGTGCAGCTGCTGTA 


CCAGACTAGCTAGAACGCCA 


60 


SSR 


Ml 04-3 


8 


104.82 


TAGAGGATCCGACGGCGT 


TGCTAGCACTCGATGAGGAA 


58 


STS 


Ml 05-3 


8 


105.12 


GTCGTCGAGTOGTCTGGAC 


TGGTCAGAGGAGAGCTAGCA 


60 


SSR 


Ml 06-1 5 


8 


106.63 


AACAAGACGGAAGCGACCA 


CATGTOTACGCCAGOTGA 


60 


SSR 


Ml 08-1 


8 


108.18 


CAATGCTGTCCGTCAATGTC 


GACATCTCGTCGATAGGCCA 


60 


SSR 


Ml 09-1 2 


8 


109.95 


GTACGTCGTGCACACCACA 


GAACGGCACCGCATGA^ 


60 


SSR 


Ml 12-5 


8 


112.05 


ATGATGTGCCTGGACCAGAG 


CCTAAGATOCOTGCTCG 


60 


SSR 


Ml 13-2 


8 


113.28 


AACACAAGCGAGGAGACGAA 


ACGATGACGACAATGGCAAG 


60 


SSR 


Ml 13-6 


8 


113.90 


ACGTCTGTCTGTGGAGTOG 


AGCAGCCTGGCAATGTOGT 


60 


SSR 


Ml 14-2 


8 


114.19 


CGTCGAGTOGCOTCATC 


GGTCACTACAAGGTCCTCGG 


60 


SSR 


Ml 14-3 


8 


114.43 


GGCATOTCGTCCTGACTGA 


TAGCACACATAGCGACATGG 


60 


SSR 


Ml 15-5 


8 


115.68 


CTAATCGTGATGGTCTCGG 


AATGACGGAATGGCAGCCTA 


60 


SSR 


Ml 17-2 


8 


117.15 


ACCTGTOATGTOCACACG 


AATGACGAGGACGGATOCC 


60 


SSR 


Ml 17-5 


8 


1 1 7.76 


GTCOTGTGOTCTGGTO 


AATCTCTCTAGCTCGTCCTCTG 


60 


SSR 


Ml 18-3 


8 


118.17 


GGTCAACGCCATCCATAACT 


CTTGCTCGTGTCGTCOTGT 


60 


SSR 


L. 1 - J 


1 




1 aA 1 L. 1 L. 1 1 1 1 UL. 


mr, A rr;rrT AT A a rrrrrrr-. 


DU 




CI -25* 


1 


34.89 


GTAGGCTCGTOGCAAAAAA 


AGAGTOAGCCGGCTATCCA 


60 


SSR 


C3-r 


3 


213.80 


CCAAGGACGCAATCTAATCG 


GTCATGGACATCGTGCTG^ 


60 


SSR 


C3-5* 


3 


214.37 


GGACAGAGCAGGTGATGTO 


GGATOGCGGACAGTOAAG 


60 


SSR 


C5-5* 


5 


17.38 


GAGGTOCACCAGTGTGCAG 


AOTCGTCCGTCCTOCTCT 


60 


SSR 


C5-24* 


5 


18.28 


GGATCGGAGGAGCCTGTOA 


TCTGTCTOTGCGTGTGTGA 


60 


SSR 


C9-2* 


9 


153.51 


TGGAGGAOTGATGTOAGG 


CTCGATGCAGTOCTOTGT 


60 


SSR 


C9-19* 


9 


153.62 


CGCAGGACATGAGGTACACC 


GCTACTCCAGTOCCAGGCA 


60 


SSR 



^Markers used for validation in BC1F2 and BC2F1 segregating populations. *Markers used for validation in a RIL population. Location corresponds to maize B73 
reference genome APGv2. 



QTL, designated qMrddl, explained 33.7% and 41.3% of 
total phenotypic variation in Jining and Taian, respectively 
(Table 4). 

Fine-mapping of qMrddl 

To fine-map qMrddl, the flanking markers umcll72 
and M103-4 were used to screen recombinants from 
BC1F2 families in Taian. Fifteen BC1F2 recombinants 
were identified and self-pollinated to generate BC1F3 
families. Of 2,685 BC1F3 individuals, 237 recombinants 
were screened and self-pollinated to generate BC1F4 pro- 
geny. To resolve recombinants associated with BC1F4 
progeny, 269 SSR primer pairs were designed within the 
qMrddl region, 34 of which were polymorphic. Finally, 
15 SSR markers (M103-4, M104-3, M105-3, M106-15, 
M108-1, M109-12, bnlgl460, M112-5, umcl858, M113-2, 
M113-6, M114-3, M115-5, M117-2, and M117-5) that 
were evenly distributed (-1-2 Mb between adjacent 
markers) throughout the qMrddl region were used to 



resolve the 237 recombinants, resulting in 23 types 
(Figure 3B, Table 2). 

Recombinant-derived BC1F4 progeny were selected 
and planted in three locations. In Taian, we grew 2,203 
BC1F4 individuals derived from 33 recombinants that 
included all 23 types. In Feicheng, we grew 2,700 indi- 
viduals derived from 37 recombinants that included 
22 types. Finally, in Jining we grew 1,805 BC1F4 indi- 
viduals derived from 31 recombinants that included 21 
types. Self-pollinated BC1F4 progeny had three geno- 
types within the heterozygous portion of the qMrddl 
locus: homozygous NT409, homozygous NT411, and 
heterozygous. DSIs for these three genotypes were separ- 
ately calculated for each BC1F4 family. For each of the 23 
recombination types, the genotype matched the pheno- 
type. Types I-VII (see Figure 3B) had the homozygous 
NT409 sequence upstream of the recombination break- 
point, and heterozygous sequences downstream. Types 
II- VII were highly susceptible to MRDD regardless of the 



Table 3 Validation of candidate regions in segregating populations 



Marker 


Location 
(Bin) 


BCi/BCiF2(Taian) 




BCi/BCiF2(Jining) 




BCi/BC2Fi(Taian) 




BCi/BC2Fi(Feicheng) 




DSI (mean 


± SD)(%) 




DSI (mean 


± SD)(%) 




DSI (mean 


± SD)(%) 




DSI (mean 


± SD)(%) 


r "VdlUc 


Homozygous 


Heterozygous 


Homozygous 


Heterozygous 


Homozygous 


Heterozygous 


Homozygous 


Heterozygous 


CI -3 


1.03 


81.19±1.03 


82.1 5 ±2.01 


0.67 


90.88 ± 1 .37 


91.19 ±0.72 


0.78 


94.02 ± 0.22 


94.1 2 ±0.25 


0.92 


76.59 ± 0.52 


75.98 ±0.73 


0.51 


CI -25 


1.03 


81.89 ±1.29 


81.06 ±1.29 


0.66 


91.09 ±0.91 


91.16 ±0.86 


0.95 


94.05 ±0.1 9 


94.08 ± 0.24 


0.94 


76.52 ± 0.60 


76.06 ± 0.52 


0.43 


C3-5 


3.08 


81.35 ±1.08 


81.87 ±1.86 


0.81 


91.26 ±0.71 


90.67 ± 1 .33 


0.69 


94.46 ±0.1 6 


93.49 ± 0.28 


0.03 


76.87 ± 0.66 


75.24 ±0.47 


0.11 


C3-1 


3.09 


81.20 ±1.04 


82.15 ± 1.90 


0.67 


91.21 ±0.70 


90.87 ± 1.41 


0.82 


94.50 ±0.1 5 


93.46 ± 0.30 


0.02 


76.69 ± 0.77 


75.67 ± 0.49 


0.20 


unncl999 


4.09 


81.91 ±1.30 


81.11 ±1.30 


0.66 


91.45 ±0.89 


90.90 ± 0.89 


0.67 


94.38 ±0.1 9 


93.76 ± 0.23 


0.14 


76.86 ±0.51 


76.70 ± 0.63 


0.81 


unncl940 


4.09 


81.91 ±1.30 


81.11 ±1.30 


0.66 


91.45 ±0.89 


90.90 ± 0.89 


0.67 


94.3 ±0.1 9 


93.78 ± 0.23 


0.17 


76.77 ± 0.62 


76.78 ± 0.52 


0.99 


umcl989 


4.09 


80.86 ± 1 .32 


82.12 ± 1.28 


0.50 


91.27 ±0.91 


91.08 ±0.88 


0.88 


94.26 ± 0.20 


93.85 ± 0.22 


0.32 


76.41 ±0.47 


76.92 ±0.61 


0.65 


C5-5 


5.03 


83.51 ±1.58 


79.81 ±1.23 


0.04 


92.54 ± 0.87 


90.11 ±0.94 


0.22 


94.59 ±0.21 


93.53 ±0.21 


0.01 


76.79 ± 0.52 


75.51 ±0.49 


0.08 


C5-24 


5.03 


83.42 ± 1 .30 


79.89 ± 1 .26 


0.06 


92.50 ± 0.93 


90.1 6 ±0.86 


0.29 


94.71 ±0.20 


93.44 ± 0.22 


2.80E-03 


76.90 ±0.51 


75.61 ± 0.68 


0.08 


bnlgl62 


8.05 


86.60 ±1.16 


77.03 ± 1 .25 


8.00E-08 


94.90 ± 0.77 


87.90 ± 0.86 


8.43E-09 


94.40 ± 0.20 


93.78 ±0.21 


0.09 


76.92 ± 0.60 


75.45 ± 0.70 


0.05 


unncl670 


8.05 


85.57 ±1.21 


77.78 ± 1.21 


1.55E-05 


93.71 ±0.84 


88.86 ± 0.88 


9.04E-05 


94.38 ± 0.30 


93.77 ±0.31 


0.16 


76.70 ±0.51 


75.64 ±0.70 


0.15 


umcll72 


8.04 


85.82 ±1.43 


77.86 ± 1 .73 


1 .04E-05 


94.07 ± 0.43 


88.45 ± 0.85 


4.98E-06 


94.42 ±0.21 


93.75 ±0.21 


0.06 


76.86 ± 0.59 


75.50 ±0.71 


0.07 


C9-2 


9.07 


83.05 ± 1 .67 


80.08 ± 1 .05 


0.11 


92.53 ± 0.88 


90.02 ±01.16 


0.06 


93.44 ± 0.26 


93.07 ± 0.25 


0.51 


76.61 ±0.68 


76.48 ±0.51 


0.89 


C9-19 


9.07 


83.05 ± 1 .67 


80.08 ± 1 .05 


0.11 


92.53 ± 0.88 


90.02 ±01.16 


0.06 


93.46 ± 0.29 


93.06 ± 0.23 


0.51 


76.61 ±0.68 


76.48 ±0.51 


0.89 
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s g g g 
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tiii 



S S 2. 2 I 2 2 

! i K 1 K ^ 
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Taian 



Feicheng 



Jining 







XIII 
XIV 
XV 
XVI 
XVII 
XVIII 
XIX 
XX 
XXI 
XXII 
XXIII 




411 



H 409 



411 H 409 



P-value No.R I 



411 



H 409 



NS 4 
NS 



NS 
NS 



NS 
NS 



S 

s 

NS 

NS 

NS 

NS 

NS 

NS 2 

NS 2 



183 46.67 

254 85.32 

92 74.07 
86 77.78 

73 74.75 
154 89.81 
60 94.07 

65 49.21 
64 45.45 

74 50.00 

93 31.11 
147 56.73 
55 52.38 
37 49.21 
58 49.21 
197 42.42 
37 41.67 
85 37.20 

66 52.05 
66 48.15 
81 51.51 
113 46.24 
63 37.37 



69.64 70.37 

83.67 82.42 

80.43 84.65 

78.47 80.95 

82.24 86.93 

85.38 86.87 

90.97 91.11 

89.50 69.70 

80.16 73.73 

83.84 88.24 

63.74 86.87 

74.60 81.48 

76.61 90.12 
82.22 77.78 
82.05 69.70 
43.02 43.79 
42.63 40.74 
37.35 38.01 
42.73 43.06 
52.08 50.33 
55.42 52.05 
47.07 43.21 
45.18 47.47 



1 .4F.-07 3 

0.54 4 

0.20 1 

0.68 NA 

0.49 1 

0.13 2 

0.32 2 

2.19E-03 1 

7.92i:-05 1 

5.05E-04 1 

6.26E-07 1 

1.02E-04 2 

6.35H-04 1 

0.02 1 

6.21E-03 1 

0.90 4 

0.87 1 

0.78 1 

0.1! 1 

0.35 2 

0.38 2 

0.32 3 



32.22 83.49 84.00 
74.30 73.74 65.71 
66.67 73.48 77.78 

NA NA NA 

88.15 82.37 82.22 

75.85 76.13 82.05 

66.67 65.97 66.01 

45.78 68.15 82.71 

36.23 64.05 77.78 
58.33 81.37 80.95 

25.49 71.85 72.22 
41.18 75.13 74.88 
33.33 57.94 63.64 
30.37 69.60 75.00 
45.78 68.15 82.72 

27.25 26.18 25.77 
40.17 40.05 54.17 
16.96 17.42 18.95 
33.33 34.51 39.26 
41.27 36.74 32.03 
33.56 33.64 38.56 

34.50 31.87 37.04 

20.26 28.37 25.93 



3.54F.-15 

0.20 

0.22 

NA 

0.31 

0.17 

0.87 
2.30E-07 
5.03i:-06 
9.16E-04 
1.59K-08 
1 .70E-08 

0.06 
2.77E-05 
6.31E-04 

0.59 

0.05 

0.91 

0.36 

0.12 

0.34 

0.32 

0.29 



4 212 

1 68 

NA NA 

1 64 

2 113 



1 66 
1 48 



2 108 
I 74 



71.85 93.78 94.44 

98.67 97.82 97.66 

98.52 96.76 96.30 
NA NA NA 

100.00 100.00 100.00 

96.14 96.07 99.03 

94.44 93.42 93.33 

71.85 96.43 93.33 

86.67 100.00 98.99 

77.78 97.85 98.41 

53.70 98.41 95.37 

62.96 96.58 94.77 

82.22 100.00 97.98 

72.84 96.33 95.56 

55.56 92.87 93.06 

63.94 63.97 65.48 

64.81 63.21 63.89 

61.90 57.08 56.94 

61.27 69.43 83.85 

86.67 79.00 79.63 

84.44 82.03 81.70 

69.59 72.86 76.00 

63.64 61.42 56.61 



i.l2i:-04 

0.47 

0.41 
NA 
NA 

0.08 

0.84 
2.27E-03 
1.19i;-04 
1.04E-04 
2.96E-08 
2.54E-12 
2.67E-05 
4.83E-06 
4.65E-06 

0.33 

0.76 

0.48 

0.06 

0.19 

0.32 

0.14 

0.18 



1.2 Mb 

Figure 3 Mapping qMrddh A. Diagram of QTL plots for qMrddl in 21 1 BC1F2 families. The logaritiim of odds profile, the relative position of 
qMrddl, and relevant markers are displayed using QTL cartographer version 2.5. B. Fine mapping of QJL-qMrddl based on recombinant progeny. 
One hundred and one BC1F3 recombinants fell into 23 types based on their genotype at 17 markers. The genetic structure of each type is shown. 
Black, gray, and white rectangles correspond to homozygous NT41 1 alleles, heterozygous NT409/NT41 1 alleles, and homozygous NT409 alleles, 
respectively. Self-pollinated progeny of these BC1F3 plants were genotyped using markers within the heterozygous region, resulting in three 
genotypes of progeny. Details concerning the number of recombinants that were planted, the number of progeny that were planted, and the 
DSI of all three genotypes of progeny are listed in the table. Significant differences {P< 0.05) among the three genotypes indicated that qMrddl 
localized to the heterozygous region and that their parental recombinant(s) was segregating (S). An insignificant difference {P >0.05) among the 
three genotypes suggested that qMrddl localized to the homozygous region and that their parental recombinant(s) was not segregating (NS). 
Analysis of both genotype and phenotype for all recombinant types enabled the positioning of qMrddl within a 1.2-Mb region between markers 
Ml 03-4 and Ml 05-3. DP: deduced phenotype, LOD: logarithm of odds. No. R: number of recombinants. No. P: number of progeny, NA: not 
available, 41 1: progeny with a homozygous genotype the same as the MRDD-resistant parental line NT41 1, H: progeny with a heterozygous 
genotype, 409: progeny with a homozygous genotype the same as the MRDD-susceptible parental line NT409. 



genotypes, whereas types I exhibited a significant differ- 
ence in MRDD resistance between the three genotypes. 
This indicated that qMrddl is located downstream of 
M103-4 and upstream of bnlgl460. Types VIII-XIII had 
heterozygous sequences upstream of the recombination 



breakpoint and homozygous NT409 sequences down- 
stream. All Types showed a significant difference in 
MRDD resistance between the three genotypes, regardless 
of the experimental location. This clearly indicated that 
qMrddl is located within the heterozygous region. Types 
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Table 4 Parameters associated with the QTL-qMrddl in the BC1F2 population 



Location 


Bins 


Flanking markers 


LOD value 


Additive 


Dominance 


(%) 


SRA 


Jining 


8.03/4 


M103-l/umcll72 


23.39 


0.14 


0 


33.72 


NT411 


Taian 


8.03/4 


M103-l/umcll72 


27.36 


0.21 


0 


41.28 


NT411 



LOD: Logarithm of odds, SRA: source of resistance allele, R^: Percentage of the phenotypic variance explained by a QTL. 



XIV and XV also showed segregation of the MRDD resist- 
ance trait and thus restricted qMrddl into the heterozy- 
gous region upstream of the breakpoint. The remaining 
types (XVI-XXIII) were resistant to MRDD regardless of 
the genotype or experimental location. This implied that 
qMrddl is located within the homozygous NT411 region 
but not in the heterozygous region. Only types XIII 
exhibited a phenotype that varied with experimental loca- 
tion. A significant difference {P < 0.05) in MRDD resist- 
ance between genotypes of types XIII BC1F4 progeny 
was detected in Taian and Jining but not in Feicheng 
(P = 0.06). This discrepancy may have resulted from 
the small number of BC1F4 progeny (43 individuals) in 
Feicheng. As such, the resistance phenotype for types 

XIII was considered to segregate, placing qMrddl 
within the heterozygous region upstream of M106-15. 
Recombination breakpoints associated with types I and 

XIV were closest to qMrddl, allowing us to fine-map 
qMrddl into the region between M103-4 and M105-3, 
a physical distance of 1.2 Mb (Figure 3B). 



Genetic model of qMrddl resistance to MRDD 

The genetic effect of qMrddl was investigated in BC1F4 
families (three genotypes and three replications). Plants 
homozygous for the NT411 qMrddl region showed a 
significantly lower DSI than the other two genotypes 
(Figure 4). In Taian, resistance to MRDD was evaluated 
for 2,203 BC1F4 plants. Plants homozygous for NT409 
qMrddl alleles had a DSI of 79.8%, which was very simi- 
lar to heterozygous plants (NT411/NT409; 77.1%). In 
contrast, plants homozygous for NT411 qMrddl alleles 
had a significantly lower DSI (49.9%). In Jining, 1,805 
BC1F4 plants were evaluated; although DSIs were gener- 
ally higher, a similar genetic pattern was observed. DSIs 
were 95.8% and 96.4% for plants carrying homozygous 
NT409 or heterozygous qMrddl alleles, respectively. In 
contrast, plants carrying homozygous NT411 qMrddl 
alleles had a DSI of 71.7%. Of 2,700 BC1F4 plants grown 
in Feicheng, plants carrying qMrddl regions that were 
homozygous for NT411, heterozygous, or homozygous 
for NT409 had DSIs of 33.2%, 67.1%, and 72.5%, respectively. 
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Figure 4 Effect of QTL-qMrddl on BC1F4 populations across three experimental sites. BC1F4 populations were divided into tliree genotypes 
(NT409/NT409, NT409/NT41 1, and NT41 1/NT41 1) according to genotypes witiiin tlie qMrddl region. Average DSI values are shown. Error bars 
indicate s.e.m. Multiple comparisons between genotypes were analyzed using SAS 9.1 PROC general linear model with Tukey's adjustment. 
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At all three test facilities there was no difference {P > 0.05) in 
DSI between heterozygous and homozygous NT409 geno- 
types at the qMrddl locus, indicating that the qMrddl QTL 
acts in a recessive manner to confer resistance to MRDD. 
Moreover, the strong genetic effect of homozygous NT411 
qMrddl alleles, which reduce DSI by 24.2-39.3%, suggested 
that this QTL represents a viable tool for enhancing maize 
resistance to MRDD. 

Validation of qMrddl in Fe RILs 

We next determined whether the qMrddl QTL is 
present in other MRDD-resistant lines of maize. Resist- 
ance was measured for 157 Fe RILs derived from the 
cross between XI 78 and HuangC. These studies were 
conducted in three locations: Taian, Feicheng, and 
Jining. SSR markers developed for the qMrddl region 
were used to screen parental lines for polymorphisms. 
Marker M103-7, which is located within the 1.2-Mb re- 
gion of qMrddl y was selected to genotype the RIL popu- 
lation. Of 157 RILs, 91 and 63 were homozygous for 
X178 and HuangC sequences at M103-7, respectively. 
Only three RILs were heterozygous at the M 103-7 locus 
and thus were excluded from subsequent analyses. RILs 
homozygous for the X178 genotype had DSIs of 30.5%, 
40.3%, and 29.3% in Taian, Jining, and Feicheng, respect- 
ively. In contrast, RILs homozygous for the HuangC 
genotype had DSIs of 56.1%, 74.3%, and 44.5% in Taian, 
Jining, and Feicheng, respectively. As such, there was a 
clear difference in MRDD resistance {P < 0.01) between 
two homozygous genotypes (Figure 5). These findings 
suggested that the qMrddl QTL could reduce DSI by 
15.2-34.0% in the RIL population, which agrees with our 
previous fine-mapping result. 




178 HuangC 

Jining 



X178 HuangC 

Taian 




X178 HuangC 

Feicheng 

Figure 5 DSI in an Fg RIL population. RILs were genotyped at 
marker Ml 03-7. DSI distributions and mean values are shown for 
two homozygous genotypes, XI 78 and HuangC at three 
experimental sites. 



Discussion 

Accurate phenotypic evaluation is critical for marker- 
trait association analyses, especially for quantitative traits 
[19]. Because large-scale inoculation of plants with 
MRDD is unfeasible and uniform infection is unreliable, 
we relied on natural infection processes. Plants were 
raised in the cities of Jining, Feicheng, Taian, and Heze 
within Shandong province, where MRDD is prevalent. 
Because of poor performance in 2011, i.e., the lack of 
planthoppers, Heze was eliminated as a testing site dur- 
ing subsequent analyses. Subsequent fine-mapping tests 
were also performed in Taian, Jining, and Feicheng to 
avoid MRDD escape. The environment significantly 
influenced MRDD development, as the disease was more 
serious in Jining compared with Taian and Feicheng. 
This may have resulted from different numbers of 
planthoppers at these three test sites. Fortunately, the 
qMrddl locus had a stable genetic effect across these 
different environments, implying that the natural infec- 
tion method was a valid approach and that our scoring 
system was appropriate for QTL analysis of MRDD 
resistance. 

Viral resistance can be influenced by both genetic 
background and the environment [20,21]. HIFs derived 
from the same cross share similar genetic backgrounds, 
making them ideal for analyzing quantitative traits. To 
identify the major QTL involved in MRDD resistance, 
50 HIFs developed from our breeding program were se- 
lected for this study. For segregating populations pre- 
pared from two HIFs with very different levels of MRDD 
resistance, a continuous distribution of resistance was 
observed rather than two distinct classes. This may have 
resulted from residual genetic-background differences or 
environmental conditions. Whole-genome SNP analysis 
revealed that 14.2% of 20,278 called SNPs differed be- 
tween susceptible (NT409) and resistant (NT411) plant 
lines. A smaller difference (3.1% of 51,628 called SNPs) 
was observed between NT401 and NT399, which are 
susceptible and resistant lines, respectively. Here we fo- 
cused on the major QTL for MRDD resistance, but 
other chromosomal regions may also be involved. A re- 
gion in chromosome 5, for example, showed a margin- 
ally significant correlation with MRDD resistance in 
both mapping populations. To generate populations for 
fine mapping we self-pollinated rather than backcrossed 
because most BC2F1 families were highly susceptible to 
MRDD, since different genotypes at qMrddl had similar 
phenotypes in BC2F1 plants in 2011. 

The recombinant-derived progeny test is an efficient 
and powerful method for fine-mapping QTLs within a 
backcrossed population when a susceptible inbred line is 
used as the recurrent parent. This method can be used 
to accurately phenotype recombinants by analyzing trait- 
marker associations in progeny [22-24]. Here we expand 



Tao et al. BMC Plant Biology 201 3, 1 3:1 45 
http://www.bionnedcentral.conn/1471-2229/13/145 



Page 10 of 13 



the application of this method to include self-pollinated 
progeny. Compared with backcrossed progeny, progeny 
generated by self-pollination can capture the effect of all 
three genotypes and create more recombinants for fine- 
mapping. However, not all recombinants produced by 
self-pollination can be used for fine-mapping, as segre- 
gating genotypes within the targeted region are not cre- 
ated when recombinants that are homozygous at both 
flanking markers are self-pollinated. Crosses that involve 
heterozygous plants represent an effective way of solving 
this problem. 

By applying the recombinant-derived progeny test to 
self-pollinated progeny, qMrddl was fine-mapped to a 
1.2-Mb region. This region decreased DSI by 24.2-39.3%. 
Forty-three additional recombinants were identified from 
8,047 BC1F5 for further fine-mapping (data not shown). 
Compared with previous reports, this represents signifi- 
cant progress towards cloning and applying qMrddl, Con- 
sistent fine-mapping results between test sites suggest that 
the recombinant-derived progeny-test represents a power- 
ful solution for fine-mapping a dominant QTL in 
backcrossed progeny or a partially dominant or reces- 
sive QTL in self-pollinated progeny. Moreover, map- 
ping results from different years (2011 and 2012) 
indicated that the genetic effect of QTL-qMrddl is 
heritable. Finally, qMrddl decreased DSI by 15.2- 
34.0% in an F^ RIL population composed of 157 fines, 
suggesting that qMrddl confers a stable genetic effect 
in diverse genetic backgrounds. 

The major QTL mapped in the current study overlaps 
with resistance regions reported by Di Renzo [17] and Shi 
[6], implying the same QTL fiinctions across different map- 
ping populations. However, this QTL has not been detected 
in the research conducted by Luan [18]. This may be due 
to a different scoring system used By Luan who adopted 
four indexes, including shorten superior internode, waxy 
enation, tassel type, and disease severity of MRDD, rather 
than an overall evaluation of MRDD symptom. 

Most important agronomic traits are quantitative in 
nature and polygenic. Compared with monogenic or 
oligogenic traits, therefore, these polygenic traits are ex- 
tremely difficult for breeders and pathologists to manage 
[21]. Isolation of QTLs, especially major QTLs, may 
simplify the analysis of quantitative traits and provide 
important resources for trait improvements. Before 
QTLs can be applied routinely to breeding programs, a 
number of chafienges must be addressed. These include 
improving diagnostic assays to detect QTLs and identify- 
ing genetic markers for marker-assisted selection [25]. In 
the present study, we used a reliable scoring system and 
developed a number of high-density markers within and 
around the qMrddl region. These tools can be used for 
widespread marker-assisted selection to improve maize 
resistance to MRDD. 



There are two categories of viral resistance in plants- 
passive resistance and positive resistance. Passive resist- 
ance is conferred by recessive plant factors, which are 
essential for the virus to complete the infection cycle. 
These typically involve protein forms that cannot be rec- 
ognized by specific viral components. In contrast, posi- 
tive resistance is conferred by dominant plant factors 
that trigger defense mechanisms in response to viral in- 
vasion [21]. The qMrddl QTL conferred MRDD resist- 
ance only in plants homozygous for NT411 alleles, 
indicating that it involves a recessive gene involved in 
passive resistance. This information wifi facilitate identi- 
fication of a candidate gene for qMrddl, 

Conclusions 

Breeding resistant maize hybrids is the most cost- 
effective way to minimize yield losses from MRDD. We 
have mapped qMrddl to a 1.2-Mb region and showed 
that it acts in a recessive manner across different genetic 
backgrounds. The discovery and fine-mapping of this 
major QTL involved in MRDD resistance lays the foun- 
dation for positional cloning of qMrddl and moves us 
closer to genetically controlling MRDD infestation dur- 
ing maize production. 

Methods 

Plant materials 

Plant materials initially selected for QTL analysis were 
50 F9 HIFs that were derived from two F5 plants of a 
single F4 individual from a hybrid, CL1165 (Figure 1). 
These 50 HIFs were evaluated for MRDD resistance in 
the summers of 2008, 2009, and 2010 at three locations, 
Taian, Feicheng, and Jining. In each location, the field 
test was conducted in randomized complete block de- 
sign (RCBD) with locations as complete blocks. 25 seeds 
from each HIF were sown in a single row 0.6 m in width 
and 5 m in length. Of the 24 F9 families with steadily 
contrasting phenotype to MRDD, two resistant (NT399 
and NT411) and two susceptible (NT401 and NT409) 
HIFs were selected to prepare two crosses, one between 
NT409 and NT411 and the other between NT399 and 
NT401. NT411 and NT409 shares 85.8% of 20,278 called 
SNPs; while NT399 and NT401 shares 96.9% of 51,628 
cafied SNPs. In 2009, two crosses were established in 
the Hainan winter nursery. In 2010, Fi plants were 
backcrossed to the susceptible parental fine in Taian. 
In the Hainan winter nursery, 211 BCi plants derived 
from the NT409/NT411 cross were self-pollinated to 
produce BC1F2 famifies. In addition, 485 BCi plants de- 
rived from the NT399/NT401 cross were backcrossed to 
NT401 to produce BC2F1 famifies. In 2011, the BC2F1 and 
BC1F2 famifies, together with parental lines, were planted 
in three locations: Heze, Taian, and Jining (or Feicheng) 
within Shandong province. In each location, the field test 
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was conducted in RCBD with locations as complete 
blocks. 25 seeds from each family were sown in a single 
row 0.6 m in width and 5 m in length. 

Based on QTL mapping, plants from the BC1F2 and 
BC1F3 families that contained recombination breakpoints 
within the target QTL region were selected for repeated 
self-pollination. In the summer of 2012, the resultant 
BC1F4 progeny were planted in three locations without 
duplication— Jining, Feicheng, and Taian. All BC1F4 plants 
were genotyped at the qMrddl locus and assayed for 
MRDD resistance. 



Survey of MRDD symptoms in the field 

Plants were grown in four locations, namely Jining, 
Feicheng, Heze, and Taian (Shandong province, China), 
and allowed to become infected with MRDD under nat- 
ural conditions, i.e., without artificial inoculations. Seeds 
were sown at May 24 to 26 to coincide with planthopper 
infestation and to make viral inoculation more likely. 
MRDD resistance was evaluated at the grain-filling stage. 
For the 50 HIFs, we scored them as resistant, intermedi- 
ate resistant, and susceptible across three years, since we 
just want to know which HIFs are consistent resistant or 
susceptible across different years at different locations. 
For QTL analysis, mapping, and fine-mapping efforts, a 
scoring system (1, 3, 5, 7, or 9) based on overall symptoms 
was adopted to evaluate MRDD resistance. Infected plants 
with different resistance scores are depicted in Figure 2A. 
The DSI was used to represent MRDD severity of families 
and was calculated as [26]: 

DSI(%) = Z (grade x number of plants in grade) 
X 100/(9 X total number of plants) 

All the phenotypic data across different replicates were 
assessed independently. 

Genotyping 

Leaf tissue was harvested for DNA extraction according 
to the SDS method [27]. SNPs were genotyped using the 
MaizeSNPSO DNA analysis kit (Illumina, San Diego, 
CA), which can survey 56,110 SNPs, using an Illumina 
BeadStation 500G at Cornell University Life Sciences 
Core Laboratories Center. Details concerning the SNP 
genotyping procedure and allele scoring have been de- 
scribed [28]. For PCR-based marker genotyping, amplicons 
were subjected to 1% agarose gel electrophoresis and visu- 
alized using a gel-imaging system (Bio-Rad Laboratories 
Inc.). Alternatively, amplicons were separated using 6% 
polyacrylamide-gel electrophoresis and visualized with 
silver-staining. 



Trait-marker association analysis 

SNPs with >50% missing data or a cluster-separation 
score of <0.3 were excluded from further analyses. TAS- 
SEL 2.0 was used to retrieve polymorphic SNPs with a 
minor- allelic frequency of >0.1, and the general linear 
model was used to analyze correlations between poly- 
morphic SNPs and phenotype. Tightly linked SNPs were 
then mapped to B73 AGPv2 [29] through BLAST com- 
parisons. Every region (<100 kb) that contained >2 co- 
segregating SNPs was considered a candidate region 
containing a QTL that conferred rice black-streaked 
dwarf virus resistance. 

Validation and mapping of the major QTL 

SSR primer pairs that covered all candidate regions were 
retrieved from public databases (http://www.maizegdb. 
org/) or developed based on B73 reference sequences as 
described by Zhang [23]. All primers were synthesized 
by Invitrogen (Beijing, China). SSR primers were first 
used to identify polymorphisms between the two paren- 
tal lines. Polymorphic SSR markers were then used to 
genotype each plant in the BCi populations. The pheno- 
type of each BCi individual was represented by DSI of 
corresponding BC1F2 families or BC2F1 families. Trait- 
marker correlations were analyzed using one-way ANOVA 
in SAS 9.1. For regions associated with MRDD resistance, 
more PCR-based markers were developed to map target 
QTLs. Linkage mapping of polymorphic SSRs was 
performed using MAPMAKER 3.0b [30]. Linkage groups 
were identified using the 'Group' command with a loga- 
rithm of odds score of >3.0. Recombination frequency was 
converted into centiMorgans using the Kosambi mapping 
function [31]. QTLs were detected using the composite 
interval mapping method [32] as with the QTL cartog- 
rapher (Version 2.5) [33]. A significance threshold for iden- 
tifying a putative QTL was obtained from 1,000 
permutations at P < 0.05 for each dataset. 

Fine-mapping of qMrddl 

The recombinant-derived progeny test [34] was used for 
QTL fine-mapping. Based on the QTL region mapped 
by winQTLcart, the BC1F2 population was screened for 
recombinants. This was followed by self-pollination in 
Shandong province in 2011. Progeny were planted in the 
Hainan winter nursery to screen for new recombinants. 
Newly-screened BC1F3 recombinants, which were het- 
erozygous at one flanking marker and homozygous at 
the other flanking marker, were selected for further self- 
pollination. This produced a segregating population for 
fine-mapping. Based on developed markers, BC1F3 
recombinants were classified into distinct types. In 2012, 
progeny of diverse BC1F3 recombinant types were 
planted in Taian, Feicheng, and Jining, with >100 kernels 
for every type in a single plot. 
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For each recombinant, the qMrddl region was sepa- 
rated into two segments, heterozygous and homozygous, 
that flanked the recombination breakpoint. Based on 
markers within heterozygous sequences, self-pollinated 
progeny were classified into three genotypes: homozy- 
gous NT409/NT409, homozygous NT411/NT411, and 
heterozygous NT409/NT411. Comparisons of score 
values between these genotypes were performed using 
one-way ANOVA in SAS 9.1. A significant [P < 0.05) or 
insignificant (P > 0.05) difference in MRDD resistance 
between these genotypes indicated that the resistance 
QTL localized to heterozygous or homozygous seg- 
ments within qMrddl y respectively. The phenotypes 
for three genotypes within the same BC1F3 recombinant- 
derived progeny were represented by DSI values. If two or 
more BC1F3 individuals shared the same donor fragment, 
they can be grouped as one recombination type. The avail- 
ability of both genotype and deduced phenotype for each 
recombinant type allowed for fine-mapping of the resist- 
ance QTL. 

Validation of the qMrddl locus in an RIL population 

The effect of qMrddl was also investigated in an Fe RIL 
population derived from a cross between XI 78 (MRDD 
resistant) and HuangC (MRDD susceptible), a commer- 
cial hybrid (ND108) widely grown in China. In 2012, 157 
Fe RILs, together with parental lines, were evaluated for 
MRDD resistance in three locations (namely Taian, 
Jining, and Feicheng) in RCBD with locations as complete 
blocks. A total of 25 seeds for each RIL were sown in a 
single row 0.6 m in width and 5 m in length. SSR markers 
generated during the fine-mapped process were used to 
genotype, i.e., screen for polymorphisms, among the 157 
RILs. Correlations between genotype and MRDD resist- 
ance were analyzed using one-way ANOVA in SAS 9.1. 
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